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Abstract 

C/3 ' We prove that the Julia set of a rational function / is computable 

Q ■ in polynomial time, assuming that the postcritical set of / does not 

! contain any critical points or parabolic periodic orbits. 

+-> . 

Lil; 1 Introduction. 

A compact subset of the complex plane is called computable if it can be vi- 
sualized on a computer screen with an arbitrarily high precision. Computer- 
| generated images of mathematical objects play an important role in establish- 

ed \ ing new results. Among such images, Julia sets of rational functions occupy 

one of the most prominent position. Recently, it was shown that for a wide 
class of rational functions their Julia sets can be computed efficiently (see 
t-h ! [6], [7], (21]) and yet some of those sets are uncomputable, and so cannot be 

| visualized (see [8], [5]). Also, there are examples of computable Julia sets 

whose computational complexity is arbitrarily high (see [3]). 
5_i ■ One of the natural open questions of computational complexity of Julia 

sets is how large is the class of rational functions (in a sense of Lebesgue 
measure on the parameter space) whose Julia set can be computed in a 
polynomial time. Informally speaking, such Julia sets are easy to simulate 
numerically. 

Conjecture. The class of rational functions of degree d ^ 2 whose Julia 
set can be computed in a polynomial time has a full measure in the space of 
parameters. 

The main result of the paper is the following. 

Main Theorem. Let f be a rational function of degree d ^ 2. Assume 



1 



that for each critical point c G Jf the uj- limit set u(c) does not contain either 
recurrent critical point or a parabolic periodic point of f . Then the Julia set 
Jf is computable in a polynomial time. 

As the reader will see below, this result can be viewed as a natural step 
towards a proof of the conjecture. 

The paper is organized as follows. In Section 1 we give all necessary 
preliminaries in Computability and Complex Dynamics, state the results of 
the paper and discuss possible generalizations. To illustrate the results of 
the paper on a simple case in section 2 we prove that for every subhyperbolic 
rational function the Julia set is computable in a polynomial time. In Section 
3 we prove the main result of the paper under a simplifying assumption that 
the rational map / does not have any parabolic periodic points. In section 4 
we complete a proof of the main result. 

Acknowledgements. It is my great pleasure to thank my supervisor 
Michael Yampolsky for posing the question of the paper and numerous fruit- 
ful discussions. 

1.1 Preliminaries on computability 

In this section we give a brief introduction to computability and complexity 
of functions and sets. The notion of computability relies on the concept of 
a Turing Machine (TM). A precise definition of a Turing Machine is quite 
technical and we do not give it here. For the definition and properties of a 
Turing Machine we refer the reader to [19] and [23] . The computational power 
of a Turing Machine is equivalent to that of a RAM computer with infinite 
memory. One can generally think about a Turing Machine as a formalized 
definition of an algorithm or a computer program. 

There exist several different definitions of computability of sets. For dis- 
cussion of different approaches to computability we refer the reader to [9]. In 
this article we use the notion of computability related to complexity of draw- 
ing pictures on a computer screen. Roughly speaking, a subset S is called 
computable in time t(n) if there is a computer program which takes time 
t(n) to decide whether to draw a given 2~™ x 2~ n square pixel in a picture 
of S on a computer screen, which is accurate up to one pixel size. Before 
giving a rigorous definition of a computable Julia set we need to introduce 
some notations. 

First we give the classical definitions of a computable function and a 
computable number. 
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Definition 1. Let S,N be countable subsets o/N. A function f : S — >■ iV is 
called computable if there exists a TM which takes x as an input and outputs 
fix). 

Note that Definition [T] can be naturally extended to functions on arbitrary 
countable sets, using a convenient identification with N. 

Definition 2. A real number a is called computable if there is a computable 
function <fi : N — > Q, such that for all n 

\a-<f>(n)\ <2~ n . 
The set of computable reals is denoted by Re- 
in other words, a is called computable if there is an algorithm which 
can approximate a with any given precision. The set Re is countable, since 
there are only countably many algorithms. The set of computable complex 
numbers is defined by Ce = Me + ^Re- Note that both Re and Cc considered 
with usual multiplication and addition form fields. Moreover, it is easy to 
see that Ce is algebraically closed. 

Let d(-,-) stand for Euclidian distance between points or sets in R 2 . Recall 
the definition of the Hausdorff distance between two sets: 

d H (S,T) = inf{r > : S C U r (T), T C U r (S)}, 

where U r (T) stands for the r- neighborhood of T: 

U r (T) = {z e R 2 : d{z,T) ^ r}. 

We call a set T a 2~ n approximation of a bounded set S if S C T and 
dH{S,T) ^ 2~ n . When we try to draw a 2~ n approximation T of a set S 
using a computer program, it is convenient to let T be a finite collection 
of disks of radius 2~ n_2 centered at points of the form (i/2 n+2 , j/2 n+2 ) for 
i,j G Z. Such T can be described using a function 

r i, if d(z,s) ^ 2~ n ~ 2 , 

h s (n,z) = l 0, if d(z, S) ^ 2 • 2~ n - 2 , (1) 
y or 1 otherwise, 

where n G N and z = {i/2 n+2 ,j/2 n+2 ), i,j G Z. 

Using this function, we define computability and computational complexity 
of a set in R 2 in the following way. 



3 



Figure 1: Values of the function hs- 

Definition 3. A bounded set S C M 2 is called computable in time t(n) if 
there is a TM, which computes values of a function h(n, •) of the form (Tjp in 
time t(n). We say that S is poly-time computable if there exists a polynomial 
p(n), such that S is computable in timep(n). 

Similarly, one can define computability and computational complexity of 
subsets of IR fc . Moreover, definition [3] naturally extends to subsets of C (see 
[H] Section 2.1). The Riemann sphere C is homeomorphic to the unit sphere 

S 2 = {x: \x\ = 1} C M 3 . 

Consider the stereographic projection 

P : S 2 \ {North Pole} C. 

The inverse of this projection is given by 

„_! f2Ke(z) 2Im(j) \z\ 2 - 1\ x 

Observe that P~ l induces the spherical metric on C, given by the formula 

dz 



Definition 4. A subset K C C is called computable in timet(n) if P 1 (K) C 
M. 3 is computable in time t{n). 
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Proposition 5. Let K C C be a bounded subset. Then K is computable as a 
subset of C if and only if it is computable as a subset ofR 2 . Similarly, K is 
poly-time computable as a subset ofC if and only if it is poly-time computable 
as a subset ofM. 2 . 

In this paper we discuss computability of the Julia sets of rational func- 
tions. For simplicity, consider the case of quadratic polynomials f c (z) = 
z 2 + c. By computing the Julia set J c of the map f c , c G C, we mean the 
following problem: 

given the parameter c compute a function h of the form (Cp) for 5* = J c . 

However, an algorithm M, computing h, can handle only a finite amount 
of information. In particular, it can not read or store the entire input c if 
c ^ Jf. Instead, it may request this input with an arbitrary high precision. 
In other words, the machine M has a command READ(m) which for any 
integer m requests the real and the imaginary part of a number <f)(m), such 
that \4>(m) — c| < 2~ m . It can be formalized using the notion of an oracle. 
Let 

D= |^:A;GZ,/gn| 

be the set of all dyadic numbers. Denote 

D" = {(xi, x 2 , . . . , x n ) : Xj G D, 1 «C j «C n} 

the n-th Cartesian power of D. We will say that the open disk Ud(c) (respec- 
tively closed disk Ud(c)) is dyadic if d and c are dyadic. Denote C the set of 
all subsets U C C such that U can be represented as a finite union of dyadic 
disks. 

Definition 6. A function <p '■ N — > D n is called an oracle for an element 
x G M. n if \\<f)(m) — x\\ < 2~ m for all m G N, where || ■ || stands for the 
Euclidian norm in R™. An oracle Turing Machine is a TM, which can 
query 0(m) for any m G N. 

The oracle (f> is not a part of the algorithm, but rather enters as a param- 
eter. In case of a computer program the role of the oracle is usually played 
by the user, who enters the parameters of the program. We should note also 
that an algorithm using an oracle <fi require m time units to read <p(rn). 

Now we are ready to define computability and computational complexity 
of the Julia set of a rational map. 
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Definition 7. Let f be a rational map. The Julia set Jf is called computable 
in time t(n) if there is a Turing Machine with an oracle for the coefficients of 
f , which computes values of a function h(n, •) of the form for S = Jf in 
time t(n). We say that Jf is poly-time computable if there exists a polynomial 
p{n), such that Jf is computable in time p{n). 

1.2 Hyperbolic maps 

A rational map / is called hyperbolic if there is a Riemannian metric \x on a 
neighborhood of the Julia set Jf in which / is strictly expanding: 

\\Df z {v)\V>\H\» 

for any z G Jf and any tangent vector v (see [H]). It follows that for a 
hyperbolic map / there is a neighborhood U of Jf on which the metric 
induced by \i is strictly expanding: 

d fl (f(x),f(y))^kd fl (x,y), for any x,yeU, where k> 1. 

Let ds be the spherical metric on the Riemann sphere. Then d^ and ds are 
equivalent on a compact neighborhood of Jf. One can deduce that there 
exists N e N, K > 1 and a neighborhood V of Jf, such that 

ds {f N (x)J N (y)) > Kd s {x,y), for any x,yeV. 

Moreover, the last property is equivalent to the definition of hyperbolicity. 
Hyperbolic maps have the following topological characterization (see [IS]). 

Proposition 8. A rational map f is hyperbolic if and only if every critical 
orbit of f either converges to an attracting (or a super- attracting) cycle, or 
is periodic. 

Braverman in [6] and Rettinger in [21] have independently proven the 
following result. 

Theorem 9. For any d ^ 2 there exists a Turing Machine with an oracle 
for the coefficients of a rational map of degree d which computes the Julia set 
of every hyperbolic rational map in polynomial time. 
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To explain why this is an important result we would like to mention the 
following. It is known that hyperbolicity is an open condition in the space of 
coefficients of rational maps of degree d ^ 2. The famous conjecture of Fatou 
conjecture states that the set of hyperbolic parameters is dense in this space. 
This conjecture is known as the Density of Hyperbolicity Conjecture. It is 
the central open question in Complex Dynamics. Lyubich [TB] and Graczyk- 
Swiatek [H] independently showed that this conjecture is true for the real 
quadratic family. Namely, the set of real parameters c, for which the map 
f(z) = z 2 + c is hyperbolic, is dense in R. 

1.3 Subhyperbolic maps 

A rational map / is called subhyperbolic if it is expanding on a neighborhood 
of the Julia set Jf in some orbifold metric (see [H]). An orbifold metric is a 
conformal metric r y{z)dz with a finite number of singularities of the following 
form. For each singularity a there exists an integer index v — v a > 1, such 
that for the branched covering z(w) = a + w u the induced metric 



in w-plane is smooth and nonsingular in a neighborhood of the origin. Douady 
and Hubbard proved the following (see |18j ) 

Proposition 10. A rational map is subhyperbolic, if and only if every critical 
orbit of f is either finite or converges to an attracting (or a super-attracting) 
cycle. 

For a subhyperbolic map, the corresponding orbifold has singularities only 
at postcritical points, which lie in Jf. In presence of these singularities the 
algorithm, which works in the case of hyperbolic maps, can not be applied 
directly for subhyperbolic maps. We show in this paper how to change this 
algorithm to compute the Julia set of a subhyperbolic map in a polynomial 
time. One of the results of this paper is the following theorem. 

Theorem 11. There is a TM with an oracle for the coefficients of a 
rational map f , such that for every subhyperbolic map f the machine 
computes Jf in polynomial time, given some finite non-uniform information 
about the orbits of critical points of f. 

We will specify later, which non-uniform information does use. 
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1.4 The main results: nonrecurrent critical orbits. 



Let / be a rational map. For a point z £ C denote 

0{z) = {f n (z) : n ^ 1} 

the forward orbit of z. Let u(z) = 0(z)\0(z) be the w-limit set of z. Denote 
Cf the set of critical points of / which lie in Jf. Put 

0/ = U w(c). 

Our main result is devoted to the class of rational maps without recurrent 
critical points. Namely, we prove here the following: 

Theorem 12. Let f be a rational map. Assume that f has no parabolic 
periodic points and Qf does not contain any critical points. Then Jf is poly- 
time computable by a TM with an oracle for the coefficients of f . 

Although in this case there is no guarantee that any kind of expansion 
holds on the whole Julia set of /, the following well-known theorem of Mane 
(see [17], [23]) implies that there is expansion on the closure of the postcritical 
set of /. 

Theorem 13. Let f be a rational map. Let M C Jf be a compact invariant 
set, such that M does not contain any critical points of f or parabolic periodic 
points and M flw(c) = for any recurrent critical point of f . Then there 
exists N EN, such that \Df N (z)\ > 1 for any z £ M. 

As we see, the theorem of Mane gives expansion only near the points, 
whose forward orbits are isolated from the parabolic periodic points. The 
question arises: is it possible to generalize Theorem [12] for the maps / with 
parabolic periodic points? It is known that the algorithms, used for comput- 
ing hyperbolic Julia sets, require exponential time in the presence of parabolic 
points (see |TS], app. H). However, in the paper [7J (see also [§]) Braverman 
proved that for any rational function /, such that every critical orbit of / 
converges either to an attracting or to a parabolic orbit, the Julia set is 
poly-time computable. Combining the algorithm used in j7] with the algo- 
rithm, which we use in the present paper to prove Theorem [T2J we obtain 
the following generalization of Theorem [T21 

Theorem 14. Let f be a rational map. Assume that Qf does not contain 
neither critical points nor parabolic periodic points. Then Jf is poly-time 
computable by a TM with an oracle for the coefficients of f . 
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1.5 Possible generalizations. 

In this subsection we discuss possible generalizations of the results of this 
paper. First we would like to mention that the algorithm which we use to 
prove Theorem [TH cannot be applied to compute the Julia set of a rational 
map /, such that Qf contains a parabolic periodic point. However, we believe 
that the following statement is true. 

Conjecture 15. Let f be a rational map. Assume that flf does not contain 
any critical points. Then Jf is poly-time computable by a TM with an oracle 
for the coefficients of f . 

Another important class of rational maps is the class of Collet-Eckmann 
maps. 

Definition 16. Let f be a rational map. Assume that there exist constants 
C, 7 > such that the following holds: for any critical point c of f whose 
forward orbit does not contain any critical points one has: 

\Df n {f{c))\ > Ce 7n for any neN. (2) 

Then we say that the map f is Collet-Eckmann (CE). The condition (TJj) is 
called the Collet-Eckmann condition. 

In [2] Avila and Moreira showed that for almost every real parameter c the 
map f c {z) = z 2 + c is either Collet-Eckmann or hyperbolic. In |T] Aspenberg 
proved that the set of Collet-Eckmann maps has positive Lebesgue measure 
in the parameter space of all rational maps of fixed degree d ^ 2. Moreover, 
there is a conjecture that for almost all rational maps / in this space the 
following is true: 

1) the forward orbit of every critical point c ^ Jf either is finite or con- 
verges to an attracting periodic orbit; 

2) for any critical point c G Jf either there exists 7, C > such that 
\Df n (f(c))\ ^ Ce in for any n e N or the forward orbit of c contains 
another critical point. 

By the theorem of Mane [131 every rational map / such that Qf does 
not contain neither critical points nor parabolic periodic points is Collet- 
Eckmann. We believe that the following generalization of Theorem [TH is 
true. 
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Conjecture 17. Let f be a rational map such that the conditions 1) and 2) 
above are satisfied. Then Jf is poly-time computable by a TM with an oracle 
for the coefficients of f . 

Thus, we conjecture that 

Conjecture 18. For almost all rational maps f of degree d ^ 2 Jf is poly- 
time computable. 



2 Poly-time comput ability for subhyperbolic 
maps. 

In this section we prove Theorem [TTJ Namely, we construct an algorithm 
A which for every subhyperbolic rational map / computes Jf in polynomial 
time. The algorithm A uses the coefficients of the map / and some non- 
uniform information which we will specify in the following subsection. 



2.1 Preparatory steps and non-uniform information. 

In this paper we actively use the classical Koebe distortion theorem (see |10j). 
Let us state it here. For 5 > 0, z G C denote 

U s (z) = {w G C : \w - z\ < 5}. 

Theorem 19. Let f : U r (a) — > C be a univalent function. Then for any 
z G U r (a) one has: 

(l-\z-a\/r)\f'(a)\ (1 + \z - a/r)\f (a)| 

{l + \z-a\/rf (i_| z _ a |/ r )3 W 

[1 + \z — a\/ry (1 — \z — a\/r) z 

The statement (HI) of the Koebe distortion theorem can be reformulated 

_ _ a\f'(a)\ 
(1+a) 2 ' '2 — (l-a)2 

U ri (a)cf(U ar (a))cU r2 (a). (5) 



the following way. Let 1 > a > 0, r\ = a }!} a ^2 , r 2 = a )!}°^2 ■ Then 



We will also use Koebe One-Quarter Theorem, which can be derived from 
Koebe Distortion Theorem. 
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Theorem 20. Suppose f : U r (z) — >■ C is a univalent function. Then the 
image f(U r (z)) contains the disk of radius \r\f'(z)\ centered at f(z). 

Fatou and Julia proved the following fundamental result. 

Theorem 21. Let f be a rational map of degree d ^ 2. Then the immediate 
basin of each attracting cycle contains at least one critical point. In partic- 
ular, the number of attracting periodic orbits is finite and do not exceed the 
number of critical points. 

In our proof of Theorem [TT] we will use the following general fact (see e.g. 
see [25]). 

Proposition 22. Let h(z) be a complex polynomial. There exists a TM 
with an oracle for the coefficients of h(z) and a natural number n as 
an input, such that outputs a finite sequence of complex dyadic numbers 
/?!, (3 2 , ■ • • , 0k for which: 

1) each f3i lies at a distance not more than 2~ n from some root of h(z); 

2) each root of h(z) lies at a distance not more than 2~ n from one of 

Denote Np(f) and Nj(f) the number of critical points of / which lie in the 
Fatou set and the Julia set of / correspondingly. The algorithm computing 
the Julia set of a subhyperbolic map will use the numbers 

AW), Nj(f) (6) 

as the non- uniform information. 

Observe that the assertions Np(f) = and Jf = C are equivalent. In- 
deed, by Proposition [TUl if Np(f) > then there is at least one attracting 
(or superattracting) periodic point and the Fatou set is nonempty. On the 
other hand if the Fatou set is nonempty, then there is at least one attracting 
(or superattracting) periodic point. The basin of this periodic point contains 
a critical point. Thus if Np(f) = 0, then the problem of computing Jf is 
trivial. Therefore, we will assume that Np(f) ^ and Jf ^ C 

Also, without loss of generality we may assume that oo ^ Jf. Indeed if 
oo G Jf, then using Proposition |22] we can find a dyadic point zq which lie 
in the Fatou set of /. Let h : C — > C be a Mobius map such that h(oo) = z . 
Consider the map h^ 1 o / o h instead of /. 

As a corollary of Proposition [22] one can obtain the following result (see 
[H], Proposition 3.3). 
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Proposition 23. Let f be a subhyperbolic rational map. There exists a 
Turing Machine with an oracle for the coefficients of f such that the 
following is true. Given the number Np(f) of critical points c ^ Jf 
outputs a dyadic set B such that 

1) all the attracting and super attracting orbits of f belong to B, 

2) for any z G B the orbit of z converges to an attracting periodic orbit, 

3) /(B)iB. 

Proof. The algorithm works as follows. Initially, let B be an empty collection 
of dyadic disks. At m-th step, m G N, do the following. By Proposition [22] 
there is an algorithm which finds all periodic points of / of period at most m 
with precision 2~ m ~ 3 . Let pi be approximate position of a periodic point of 
period fcj. The corresponding periodic point z% belongs to the disk U 2 -m-3(pi). 
Consider the disk U 2 - m /2(pi). If U 2 -™./i{pi) does not intersect neither one of 
the disks from B, then approximate the image f ki (U 2 - m /2(pi)) by a dyadic 
set with precision 2 _m_3 . Namely, find a set Wi G C such that 

d H (m,/*(tf 2 -v a (Pi))) <2- m ~\ 

Verify if 

[/^(^CC/^fe). (7) 

This would imply that 

f h {U 2 - m/ 2{pi)) C U 2 - m -l{Wi) C U 2 -m/2_ 2 -m-l(pi). (8) 

In this case, compute dyadic sets Bj, such that B = U 2 -m(Wi) and for each 
j = 0, 1, . . . , hi — 1, f(Bj) d Bj + i, where j + 1 is taken modulo fcj. Add 
dyadic sets Bj to the collection B. 

Next, calculate approximations of the images f m (ci) of critical points 
of / such that |sj — f m (ci)\ < 2 -m ~ 1 . If there are N F (f) of the points Sj, 
such that 

U 2 -m-l( Si ) CB=US, 

then we stop the algorithm and output B. Otherwise, go to step m + 1. 

Let us show that the algorithm eventually stops and outputs a set B, 
satisfying the conditions 1) — 3) of Proposition |22J Let z be an attracting 
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(or super-attracting) periodic point of period k with multiplier A. Let |A| < 
r < 1. Then for small enough e > one has 

f ki (U £ ( Zi )) c U rE { Zi ). 

It follows that for some m > k the corresponding approximation pi of z and 
the set Wi satisfy the property ((?]). On the other hand if ([8]) holds, then 
by Schwartz Lemma U 2 -m/2(pi) contains an attracting periodic point, whose 
basin contains U 2 -m/2(pi). Therefore if the algorithm runs sufficient amount 
of steps, the union B of dyadic sets from B satisfies the condition 1) of the 
Proposition [231 

By Proposition [10J the orbit of each critical point of / which does not lie 
in Jf converges to an attracting periodic orbit. Thus, for some m there will 
be Np(f) of the points Sj, belonging to B. This implies, then the algorithm 
stops. By Theorem [2H we obtain that B contains all attracting periodic 
orbits of /. Notice that conditions 2) and 3) of the Proposition I2U1 for this 
set B are satisfied by construction. □ 

Denote CFf the set of critical points of / which lie in the Fatou set of /. 
Put 

PF f = \Jf j (CF f ). 

j>o 

The next statement is a subhyperbolic analog of Proposition 3.7 from [9]. 

Proposition 24. There exists an algorithm which, given the coefficients of a 
subhyperbolic rational map f of degree d ^ 2 and the number Np{f), outputs 
a planar domain U G C such that: 

(1) U € f{U), 

(2) f{U) n PF f = 0, 

(3) J f m U. 

Proof. First use the algorithm from Proposition (T23]) to find a dyadic set B 
satisfying to the conditions 1) — 3) of Proposition [221 Let m G N. We can 
algorithmically construct a dyadic set W such that 

f- m {B) m W m / 1_m (5). 
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Compute a dyadic number d > such that 

U d {W)mr\W). 

Also, compute dyadic approximations Sj of critical points c, of / such that 

| Sj — Ci\ ^ d. 

By Theorem [TQ], for any critical point q G C-F 1 / we will eventually have: 

Therefore, for large enough m, the set W will contain Np(f) of the points 
Sj. Take such m. Then for any S, G one has 

Ci G G t^(W) C r 1 ^). 

Thus, CFf C Compute a dyadic set W such that 

f~ 3 (W) CW C r 2 (W). 
Clearly, for the set U = C\W conditions (1) - (3) hold. □ 

2.2 Construction of the subhyperbolic metric. 

Here we give the construction of a subhyperbolic metric (see [H]). We modify 
the construction from [18J to be able to write it as an algorithm. First we 
recall the definition and basic properties of an orbifold. We refer the reader 
to [IB] for details. 

Definition 25. An orbifold (S, v) is a Riemann surface together with a func- 
tion v : S — > N such that the set {z G S : v(z) ^ 1} is discrete. Points z for 
which u(z) 7^ 1 are called branch points. 

Definition 26. Let f be a subhyperbolic rational map. An orbifold metric 
fi on a neighborhood U of Jf is called subhyperbolic if f is strictly expanding 
on U with respect to fi: 

PV(*)IU>A>1 
for any z G f^iJJ) except the branch points. 
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Let p : S' — > S be a regular branched covering. Then for every z G S 
the local degree of p at a point w G p -1 ^) does not depend on w. One can 
define the weight function v : S — > N of the covering p assigning z/(z) the 
local degree of p at to G f" 1 {z). 

Definition 27. Lei (S", v) be an orbifold. A regular branched covering 



with the weight function v such that S' is simply connected is called a uni- 
versal covering of the orbifold (S, v). We will use the notation S v — > (S, v) 
for a universal covering of this orbifold. 

Proposition 28. Let (S, v) be an orbifold. The universal covering 



exists and unique up to conformal isomorphism, except in the following two 
cases: 

1) S » C (the Riemann sphere) and S has only one branch point; 

2) S ~ C and S has two branch points a±, 02 such that v{a\) 7^ 2/(02) • 

The Euler characteristic of an orbifold (S, v) is the number 



Since the set of branch points of S is discrete the last sum contains at most 
countable number of nonzero terms. If S contains infinitely many branch 
points then we set x(S, v) = —00. The orbifold (S, v) is called hyperbolic if 



Lemma 29. If (S, v) is a hyperbolic orbifold then S u conformally isomorphic 
to the unit disk. 

Let / be a subhyperbolic rational map. Let U be the set from Proposition 
EU Construct an orbifold (U, v) in the following way. Put S = U . Denote 
CJf the set of the critical points of / which lie in Jf. As the set of branch 
points of U take BP = {p(c),c G CJf,j G N}. Since / is subhyperbolic, 
BP is finite. Put v(z) = 1 for all z G U\BP. Denote n(f, z) the local degree 



p : S' S 




X (S,v)<0. 
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Figure 2: Construction of the subhyperbolic metric 



of / at z. Define numbers u(a),a G BP such that the following condition 
holds: 

for any zEU v(f(z)) is a multiple of v(z)n(f,z). (9) 

If the orbifold (U, v) that we obtained is not hyperbolic, take any repelling 
orbit {zj} of / of the length at least 5 and replace f(zj) with 2v(zj) for all 
Zj from this orbit. The new orbifold will be hyperbolic and satisfying the 
condition (jUj). 

By Proposition |5S] there exists a universal covering 

U u ^(U,u). 

By Lemma [29j without loss of generality we may assume that U u = U is the 
open unit disc. Since V = f~ l (U) C U, Condition (jHJ) guarantee that the 
map f~ l lifts to a holomorphic map 

F : U ->■ U. 

Note that W = F(U) = 7r _1 (\/) is strictly contained in U. Denote pu the 
Poincare metric on U. By the Schwartz-Pick Theorem, the map F is strictly 
decreasing in the metric pu- Let \i be the projection of pu onto U. For any 
ioGIJ and z = tt(F(w)) one has: 

\\Df(z)\\,= \\DF(w)\\-\ (10) 
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It follows that the map / is strictly expanding with respect to the norm in- 
duced by /i. To show that the metric /i is subhyperbolic we need to prove that 
jj, is uniformly strictly expanding. First we will prove an auxiliary lemma. 

Lemma 30. There exist a constant C > such that for any z G U except 
the branch points of the orbifold (U, v) one has: 



< 



fi{z) 



dz 



< C max 



where aj are the branch points. The constant C can be obtained construc- 
tively. 

Proof. Let z G U. Assume that z is not a branch point. To estimate /j,(z) 
without loss of generality we may assume that z = vr(0). Recall that the 
Poincare metric on the unit disk is of the form: 



Therefore one has: 



2\dw\ 



\w\ 



dz 



Dtt 



We can construct a dyadic number R such that U is contained in a disk of 
radius R. By Schwartz Lemma, 



Thus, 



|L>7r(0)| ^ R. 
fi{z) 



dz 



> 2R- 



Recall that the set BP of branch points consists of critical points which lie 
in Jf and possibly one repelling periodic orbit. It follows from Proposition 
22] that we can construct a dyadic number e > such that disks C/ 2£ (cij) 
are pairwise disjoint and all belong to U. Assume first that z ^ [jU e (aj). 
Considering the branch 

6 : UJz) ->U, 



of 7r 1 such that <f>(z) 



0, by Schwartz Lemma we obtain 
\Dn(0)\ > e. 
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\ 




h(w)=w m +aj 



Figure 3: Illustration to Lemma [301 



Let z G U e {dj) for some j. Then in a neighborhood V(0) of the map 7r can 
be written in the form 

ir{w) = g{w) m + a, h 

where giw) is a one to one map from V^(0) onto Us(0) with 5 = (2s) l l m . In 
V(0) one has 

Dn(w) = mg{w) m - 1 Dg{w). (11) 
Consider the map x — 9 1 '■ U$(0) — > V(0). Since 

|^(0)| = \z- aj \ 1/m < 2~ 1/m 5, 
by Schwartz Lemma we get: 

\D X (g(0))\<K = K{e), 
where constant K > can be obtained constructively. Now (fTTl) implies 

|£>tt(0)| ^ mir 1 |2-a i | 1 - 1/w , 
which finishes the proof. □ 
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Proposition 31. There exists a constant A > 1 such that 

P/(z)IU>A 

for any z G V = f (U). The constant A can be constructed algorithmically. 

Proof. For a map g : U\ — > U2 between two hyperbolic Riemann surfaces 
denote 

\\Dg(z)\\ UuU2 

the magnitude of the derivative of g computed with respect to the two 
Poincare metrics. Denote dist M the distance in the metric \i on U and distu 
the distance in the Poincare metric on U. Let 

% : W 

be the inclusion map. Then we have: 

\\DF{w)\\ u,u — ll-D-f 1 ( u ')l|u,vy||-D^(-P 1 (w))||Ty,u ^ \\Di{F(w))\\w,v- 

Using Lemma [30] it is not hard to show that we can construct a dyadic 
constant R > such that 

dist^O, U \ V) < R for any z e V. 

Let z G V . Then there exists w G W and ( G U \ W such that 

7r(u)) = z and distu(a;, Q < R. 

A suitable fractional linear transformation sends ( to and u to x > 0. 
Explicit calculations show that 

1 + x e d — 1 e R — 1 

d = d v (x,0) = log- , so that g = - ^ - . (12) 

1 — x e a + 1 e K + l 

Now, by Schwartz-Pick Theorem 

||-D*Mllw,u ^ \\Di(x)\\ v \ {0}!V . (13) 

The right hand side of the last inequality can be estimated explicitly. It is 
equal to 

2|xlogx| 
a ^ = !_ X 2 < L 

Note that a(x) increases with x. By f|T2l and ([TBI the value 1 1 (if) 1 1 w,u is 
bounded from above with a(X) < 1 for X = (e R - l)/(e R +1). By ([TO]) we 
obtain that ^/(z)^ for z G V is bounded from below with l/a(X). □ 
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Corollary 32. The metric fi is subhyperbolic. 

Lemma EDI and Proposition EH together make a subhyperbolic analog of 
Proposition 3.6 from [9]. 

2.3 The algorithm. 

Denote V k = f~ k {V). Notice that 

J f 14+i mV k mV 
for any fceN. Let us prove the following auxiliary statement. 

Proposition 33. There is an algorithm computing two dyadic constants 
Ki,K 2 >0 such that for any z G V 3 \ Jf and any k 6 N if f k (z) G Vx \ V 3 
then one has 

Kx < n . K 2 
\DP(z)\^ d[Z ' Jf) ^\Df»(z)\- 

Proof. First construct a dyadic number R such that 

< R < min{rf(C \ V 3 , J f ), d(V 1} C \ V)}. 

Then for any z, k, satisfying to the conditions of the proposition, the open 
disk UR(f k (z)) does not intersect neither Jf nor a forward orbit of a critical 
point of /. Let 

: U R {f k {z)) -> C 

be the branch of {f k )~~ l such that <p(f k (z)) = z. Then by Koebe Quarter 
Theorem the image 4>(Uit{f k {z))) contains the disk or radius ^R\D<j)(f k (z))\. 
Since Jf is invariant under / it follows that 

Set Kx = R/4. 
Further, denote 

PJ f = {P(c):ceCJf,j>0}. (14) 

Notice that PJf is finite. Recall that we can algorithmically construct ap- 
proximate positions of all critical points of / which lie in Jf with any given 
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Figure 4: Illustration to Proposition l33l 



precision. Thus, we can approximate PJf- We can also algorithmically con- 
struct positions of some points in Jf which lie apart from PJf with any 
given precision. For instance, using the algorithm from Proposition [22] we 
can calculate approximate position of a repelling periodic orbit. Using the 
above we can construct a finite number of pairs of simply connected dyadic 
sets Wk <s £4 such that the following is true 

1) Wj H J f ^ for any j; 

2) \jWjDVi\Vs; 

3) \jUjCV\PJf. 

For each j fix the Riemann mapping ipj : Uj — > D = {z : \z\ < 1}. Assume 
now that z, k satisfy the conditions of the proposition. Then f (z) C Wj for 
some j. Let : Uj —> U(z) be the branch of (f h )~ 1 such that cj)(f k (z)) = z. 
Consider the map 

(f) o ijjj 1 : D — > C/(^). 
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Notice that ipj(Wj) <s D. Applying both parts of Theorem [T§1 to the map 
o ipj 1 , we can construct a dyadic number r 3 - > not depending either on 
z or on fc such that </>(Wj) is contained in the disk of radius rj\D<ft(f k (z))\ 
centered at z. It follows that 



d(z,J f )< '' 



\Df k {z)Y 

Set K 2 = max{rj}. □ 
As a preparatory step construct a dyadic set W 2 such that 

V 3 mW 2 V 2 . 

Compute dyadic numbers s, e > such that 

e < mm{d(W 2 , C\V 2 ), d(V 3 , C \ W 2 )}, V 1 C U s (J f ). (15) 

Let m = max{i/(aj)}. Lemma [301 implies that for any z G V 

C~ x d{z, Jf) < dist M (z, J/) < mCd(z, Jf) 1/m . (16) 

Let log stand for the logarithm with base 2. We will use the standard notation 
[x] for the integer part of a real number x. 

Assume that we would like to verify that a dyadic point z is 2~™~ 1 close 
to Jf. Consider first points z which lie outside V 3 . Construct a dyadic set 
W 3 such that 

Jf C Wz <& V 3 . 

Then we can approximate the distance from a point z £ W3 to Jf by the 
distance form z to W$ up to a constant factor. 

Now assume 2 e V3. Consider the following subprogram: 
i := 1 

while i ^ [log(mCV /m )/ log A + (n + 1)/ log A] + 1 do 
(1) Compute dyadic approximations 

Pi « /*(*) = /(Z* -1 ^)) and d, » |£>/*(z)| = iD/^z) • Df{f-\z))\ 



with precision min{2 n 1 ,s}. 
(2) Check the inclusion G W 2 : 
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• if pi eW 2 , go to step (5); 

• if pi W 2 , proceed to step (3); 

(3) Check the inequality di ^ K 2 2 n+l + 1. If true, output and exit the 
subprogram, otherwise 

(4) output 1 and exit subprogram. 

(5) i -> i + 1 
end while 

(6) Output end exit, 
end 

The subprogram runs for at most L = [log(mC 2 s 1 /' m ) / log A+(n+l) / log A] + 
1 = 0(n) number of while-cycles each of which consist of a constant number 
of arithmetic operations with precision 0(n) dyadic bits. Hence the run- 
ning time of the subprogram can be bounded by 0(n 2 lognloglogn) using 
efficient multiplication. 

Proposition 34. Let f(n,z) be the output of the subprogram. Then 
( 1, if d(z,J f ) >2- n ~\ 

I 0, if d(z,J f ) <K2- n ~\ (17) 

y either or 1, otherwise, 

where K = : , 

Proof. Suppose first that the subprogram runs the while-cycle L times and 
exits at the step (6). This means that pi G W 2 for i = 1, . . . , L. In particular, 
p L £W 2 . It follows that f L (z) £ V\. By $53} and we obtain: 

d(z, J f ) <: Cdist^, Jf) ^ C\- L dist^f L (z), J f ) ^ 
\- L mC 2 d(f L (z), J f ) 1/m ^ X~ L mC 2 s 1/m ^ 2~ n ~ 1 . 

Thus if d(z, Jf) > 2~ n ~ 1 , then the subprogram exits at a step other than (6). 

Now assume that for some i ^ L the subprogram falls into the step (3). 
Then 

Pi-i e W 2 and pi W 2 . 

By (J35J, f\z) G Vi \ V 3 . Now if d, t ^ K 2 2 n+1 + 1, then \Df(z)\ ^ K 2 2 n+1 . 
By Proposition [331 

d{z, J f ) ^ 2- n -\ 
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Otherwise, \Df(z)\ < K 2 2 n+1 + 2 ^ (K 2 + l)2 n+1 . In this case Proposition 
implies 



d{z,J f )>—^2 



□ 



Now, to distinguish the case when d(z, Jf) < 2~ n ~ 1 from the case when 
d(z, Jf) > 2~ n we can partition each pixel of size 2~ n x 2~ n into pixels of size 
[2r n / K) x (2 _n /K) and run the subprogram for the center of each subpixel. 
This would increase the running time at most by a constant factor. 

3 Maps without recurrent critical orbits and 
parabolic periodic points. 

In this section we will prove Theorem[T2J Throughout this section let / stand 
for a rational map without parabolic periodic points such that flf does not 
intersect the set of critical points of /. 

3.1 Preparatory steps and nonuniform information. 

As in the case of subhyperbolic map, without loss of generality we will assume 
that oo ^ Jf. For any c G CJf put 

iVo(c) = max{n : f n (c) eCJ f } + l. 

Put iVo = max iVo(c). Denote 

ceCJf 



C = {f n (c) : c e CJ f , 0^n< N (c)}, M = {f n (c) : c G CJ f , n ^ N (c)}. 

By our assumptions, there are no either recurrent critical orbits or parabolic 
periodic points of /, the set M is invariant and does not contain critical 
points of /. Thus, the set M satisfies the conditions of Mane's Theorem 
f fT3|) . The following result is classical (see [IS]). 

Theorem 35. Let g be a rational map. Then the boundary of each cycle of 
Siegel disks and each cycle of Herman rings belongs to PJ g (see (fl^PJ. 

Lemma 36. There are no either Siegel disk cycles or Herman ring cycles in 
the Fatou set of f . 
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Proof. Assume for simplicity that there is a Siegel disk A. By replacing 
/ with an iterate if necessary, we can assume that /(A) = A. Then the 
boundary dA of the Siegel disk is forward invariant under / and belongs 
to PJf. It follows that dA does not contain any critical point of /. Thus, 
dA satisfies the condition of Mane's Theorem ffl3l) . Therefore, there exists 
N such that f N is expanding on a neighborhood of dA. This is impossible 
since / is conjugated to a rotation inside A. The other case can be treated 
similarly. □ 

To compute the Julia set the algorithm will use the following non-uniform 
information: 

NT. N F (f), Nj(f) and degrees mi,...,mjv J .(/) of the critical points of / 
which lie in Jf, 

N2. N , JVeN, dyadic numbers 5, 5' > 0, q > 1 and a dyadic set U D M 
such that 

U S/ 2(M) Df N (U), UdU s ,(M), U 5 (U)nC = 
and for any z G Us(U) one has 

\Df{z)\>q. 
In this section we will prove the following theorem. 

Theorem 37. Let f be a rational map such that f has no parabolic periodic 
points and Qf does not contain any critical points. There exists a Turing 
Machine which, given an oracle for the coefficients of the map f and the 
non-uniform information (N1,N2), computes Jf in a polynomial time. 

Now we prove several auxiliary lemmas. 

Lemma 38. For any z G U and any z±,Z2 G Us(z) one has 

\f N ( Zl )-f N (z 2 )\ >q\ Zl -z 2 \. 

Notice that, in particular, the restriction of / on Ug(z) is one to one for 
any z G U . 

Proof. Using Lagrange formula 

g(zi) - 9(z 2 ) = (*i - z 2 )Dg(\ Zl + (1 - \)z 2 ), A G [0, 1], (18) 
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h 

Figure 5: Illustration to Lemma [391 

for g(z) = f N {z) we obtain 

\f( Zl )-f( Z2 )\ >q \ Zl - Z2 \. 

□ 

For simplicity set F = f . Lemma [381 implies that 

d{F{z), J f ) ^ qd{z, J f ) for each z G U S {U). (19) 

Lemma 39. For any z G U and k ^ min{j G N : F^(z) ^ U}, k G N, one 

has 



l -\DF k {z)\d{z,J f )^d{F k {z),J f ) ^ 9 -\DF k (z)\d(z,J f ), 

d{F k {z),M) <: ^ \DF k (z)\ d{z,M). 

Proof. Let z, k satisfy the conditions of Lemma |39l Put 

z j = F j {z),U j = U s {z j ),j = 0,1,..., k. 
By Lemma [Ml F is univalent on Uj for each j = 0,1,..., k — 1 and 

F(Uj) D 

It follows that there exists V C Uq, such that 

F k (V) = U k 
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(20) 



and F k is univalent on V. Lemma [2H1 implies that 



V c U s/q (z ). 

Let h : U k —> V such that ho F k is identical on V. Put 



r = d(z k ,J f ). 



Then r ^ |. 



By Koebe Distortion Theorem (see (jSD), one has: 



^gr|fc'(* fc )|( Z ) C h i U r{Zk)) C £/ 2 r|fc'(* fc )|0)- 

Since J/ is invariant under /, it follows that 

\r\h\z k )\ ^ d(z,J f ) ^ 2r\h'(z k )\. 



(21) 



(22) 



Since Dh(z k ) = D Jk^ , ( |22|) is equivalent to the first part of ( 120]) . Similarly, 
one can prove the second part of ( 1201) . Observe that we can not prove a lower 
bound for d(F k (z), M) in this way since M is not necessarily invariant under 



Till the end of this section in each of the preceding statements we as- 
sume that the Turing Machine has an access to the non-uniform information 
(Nl, N2) and an oracle for the coefficients of the map /. 

Lemma 40. There exists a TM which computes dyadic numbers e\ > 0, K > 
0, such that for any critical point c of f which lie in Jf one has 



for any z 6 U ei (c), where m is the degree of f at c. 

Proof. By Proposition [22] we can find approximate positions of all critical 
points of /. Notice that the orbit of every critical point of / which lie in the 
Fatou set converges to an attracting cycle. Using the ideas from the proof 
of Proposition [23] we can distinguish the critical points which belong to the 
Fatou set from the critical points which belong to the Julia set. Let c G Jf 
be a critical point of degree m. Then we can approximate the coefficient a m 
at (z — c) m in the Taylor expansion of / at c and compute 7 > 0, such that 



F-\ 



□ 



d(f(z),J f )^Kd(z,J f )\z-c 



f(z) - f(c) 
a m (z - c) m 



- 1 < 1 
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for any z G U 1 (c). Let oo be any m-th root of a m . Then there exists a unique 
holomorphic map ip from £/ 7 (c) to a neighborhood of the origin, such that 

f( z ) = ^( z ) m + /(c) on C/ 7 (c) and ^'(c) = w. 

Make 7 smaller if necessary in order to have 

\i//(z)\ ^ \u\/2 for any z G C/ 7 (c). 

The image /(C/ 7 (c)) contains a disk U a (f(c)), where a > can be algorith- 
mically constructed. Let < S\ < 7, such that 

f(U £l (c))GU a/2 (f(c)). 

Take 2 G U ei (c). Observe that 

d(f(z),J f )^d(f(z),f(c))^a/2. 

Let (0 E Jf, such that 

d(/(*),J/) = |/(*)-Co|. 

Then ( G U a (f(c)). There exists z G J/ fl C/ 7 (c), such that /(zo) = Co- One 
has 

m— 1 

/(*) - f(z ) = ^(z) m - ^(z ) m = n - e 2 ^/-V(^o)) • 

3=0 

Let ^ jo ^ m — 1 be the number for which e 2mjo ^ m "ip(z ) takes the closest 
value to ip(z). From simple geometric observations it follows that 



Wz) - e 2 ^ /m i){z Q )\ > Wz)\ sin(7r/m) > ^ M 

1 1 mm 



z — c 



for any j 7^ j - Since *0 1 (e 27 ™- 70 /™"^'^)) G J/, it follows that 

(I II I \ TO— 1 

^ ) W*) - e 2 ^/^(^ )| > K\z - c\ m - l d(z, J f ), 

where K = m 1 ' m \uj\ m /2. □ 
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Using the definition of the dyadic set U and Proposition [22] one can prove 
the following statement. 

Lemma 41. There exists a TM which computes a dyadic number 62 > 
such that 



u £2 (c) n f N (u) = and f No (u E2 {cyj C u. 



(23) 



Proposition 42. There exists a TM which computes dyadic number < 
7 < min{ei, 62] such that the following is true. For any z £ U 7 (^C^j and 
k = min{j £ N : f Nj+N °(z) £ U} one has 

d(f kN+N °(z),J f )>qd(z,J f ). 

Proof. First put 7 = min{ei,£ 2 }- Let z £ U 1 (Cj \ Jf. Set 

w = f No (z), k = mm{j £ N : F j (w) £ U}, 

where F = f . Using Lemma |39| we obtain: 

d(F k (w),J f ) > d(F k (w),M) 

d(w,J f ) ' 9d(w,M) ' 1 ' 

Let c £ C. Put h = f N °. Let m c be the degree of h at c £ C. Then we can 
algorithmically construct a dyadic number L c > 0, such that 

\h(z)-h(c)\ < L c \z-c\ mc 

for any 2 £ i7 7 (c). Then 

M ) < <*(M*)> ^( c )) < ^l- 2 - c \ mc - ( 25 ) 

Lemma HQ] implies that we can algorithmically construct dyadic numbers 
a c > 0, K c > such that 

d(h(z), J f ) ^ K c d{z, J f )\z - c\ m '- 1 (26) 

for any z £ U ac (c). Combining (125'j) with (l2B"j) . we get 

— — — ^ K C L C — : — (27) 

a(z,Jf) \z — c\ 
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for any z G U min{ac ^ } (c) \J f . 

Take 7 such that 7 min | - K ^ c , a c j for any c E C. Recall that in 
(124)) F stands for and w stands for /^"(s), in ^ h stands for J^ . Now 



( 1271) implies that 

d{z, Jf) 0' 

for any z G ?7 7 (cj- Combining ( |24|) and (133|) . taking into account that 
U D Us'(M), we obtain the statement of Lemma 1421 □ 

A direct analog of the Proposition 1241 holds for the maps /, which satisfy 
the conditions of Theorem [T2j 

Proposition 43. Let f be a rational map such that f does not have any 
parabolic periodic points and Qf does not contain any critical points of f. 
There exists a TM which given an oracle for the coefficients of the map f 
and the non-uniform information outputs a planar domain V G C such that: 

(1) V m f(V), 

(2) p{V) n PF f = 0, 

(3) Jf m V. 

Denote 

& = CJ f U |J 0(c) = CUM, W = V\&. 
cecjf 

If CJf = then / is hyperbolic. Since for hyperbolic maps poly-time com- 
putability of the Julia set is well known (see [6] and [21]), we will assume 
that CJf 7^ 0. Then /(O) is strictly smaller than G. It follows that f{W) is 
strictly larger than W . Let || • ||w stands for the hyperbolic norm associated 
with W . Then 

\\Df(z)\\ w > 1 for every z G W. 
Let d w be the metric on W induced by || • \\w It follows that 

dw(f(z), Jf) > d w (z, J f ) for every z G f~\W) n W. (29) 
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Notice that 

inf \f'(z)\>0. 
*eJ/\(Oy(o)) 

One can construct a dyadic subset V\, such that 

r\V) m V 1 m V. 

Denote W 1 =V 1 \G. 

Lemma 44. One can algorithmically construct dyadic numbers r G (0, 1], t > 
1, -R > 0, e > 0, smc/i that 

(1) d w (f(z), J f ) > td w (z, J f ) for anyzeW\ (V 7 (<5) U tf); 

^ Jf) > rd(z, J f ) for any zeW\U^ ; 

(3j D R (0) DWDWrD U e (J f ). 

Notice that the Euclidian metric d and the hyperbolic metric dw are 
equivalent on any compact subset of W. The following lemma can be proven 
similarly to Lemma 1301 

Lemma 45. One can algorithmically construct a dyadic constant C > such 
that 

C- X d{z, Jf) ^ d w (z, J f ) ^ Cd{z, Jf) (30) 
for any z G W\ \ (\J 1 (d^j U if) . 

Lemma 46. There exists a TM which computes dyadic numbers S, K > 
such that for any k G N and z G C one has d (f k (z), JA ^ Sd(z, J f) K ■ 

Proof. Since / does not have parabolic periodic orbits one has 

mfd(P(C\V l ),Jf)>0. 

Therefore there exist S < 1 , Kq > 1 such that for any z £ V% and k G N 

d(/ fe (z),J/) >Sd{z,J f ) Ktt . 
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Make S smaller if necessary in order to have 



d(C \ V u J } ) > S sup d(z, J f ) Ko . 



z£Vi 

Let zq G Vi \ Jf. Assume there exists / G N such that 

d(f(z Q ),J f )^Sd(z ,J f ) K °. 
Let I be the minimal such number. Then 

P(z )eV 1 foij = o,i,...,L 

It follows from Lemma HD] that one can compute dyadic numbers 1 ^ 77 > 0, 
m ^ 1 such that 



d{f j (z),J f ) >vd(z,J f Y 



(31) 



for any z G U 1 ICj and j = 1, 2, . . . , N + N - 1. Put Z = 0. Denote 
inductively /j and Zi = f li (zo) as follows: 



if m & [C U 17. 



— < 



/i + TV, if ^ G C/ and li + N ^ I, 

k + JV + fciV, if Zj G U 7 (c) and Z< + 7V + fciV ^ /, 

L otherwise. 



(32) 



while U < I, where 

k = min{j G N : f Nj+N °(Zi) $ U}. 
Now by (I29p . Lemma HH and Proposition 1421 the following is true: 

1) d w (z i+1 ,J f ) ^ J/) for any i; 

2) if zj ^ U y (c) U £/, then 

d\v(zi+i,Jf) ^ td w (z u J f ) and d(^+i, J/) ^ rd(zj, J/); 

3) if Zi G £/ and k + N ^l, then 

d(z i+1 ,J f ) ^ gd(zj, J/); 
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4) if Zi e U and k + N > I, then 

d(z i+1 ,J f ) ^ r N ~ 1 d(zi, Jf); 

5) if Zi e U 1 (dj and k + N + kN < /, then 

d(z i+1 ,J f ) ^ d(zi,J f ); 

6) if Zj G I7 7 (cj and h + N + kN > I, then 

d(z i+ i, Jf) ^ r]r N ' 1 d(z i , J f ) m . 

The items 1) — 5) are direct corollaries of the choice of the dyadic constants, 
formula fl29|) . Lemma HH and Proposition H2J Let us proof 6). Let 

Zi e C/ 7 (dj and h + N + kN > I. 

If ij + No > I, then 6) follows directly from (I3TI) . Otherwise denote n the 
remainder of I — Zj — N modulo N. Then by (1311) and (TT9|) 

d(z i+1 ,J f ) >r N - l d(f l - n (zi),J f ) >r N ~ x d(f li+N \z i ),J f ) ^r^V^^/)" 
Let s be the number of i-th, for which 

Zi £ U 1 (c\ U U. 

If s 7^ then let Ji be the minimal and I2 be the maximal such i. 
case 1: s = 0. Then one of the following two possibilities holds. 

a) z E U 1 (d\ Then I = h. By 6), 

d(f l (z ),J f )>r,r N - l d(z ,J f ) m . 

b) z e U. Then by 3) and 4), 

d(f\z ),J f )^r N - 1 d(z ,J f ). 



case 2: s ^ 0. There are two possibilities, 
a) Z/ 1 = /. Then 7j = J 2 - As in the case 1, 

d(f l (z ),J f )> V r N - l d(zo,J f 
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b) li x < I. Then, by 3) and 5), d(zi 1 , Jf) ^ d(z , Jf). By (13"U1) . 

d(z l2 , J f ) > C- l d w (z l2 , J f ) > C~ x d w {z Ix , Jf) > C- 2 d(z , J f ). (33) 

Among the numbers 1% + 1, 12 + 2, . . . , I there is at most one number j such 
that 

Zj e U 1 (c 
and at most one number j such that 

Zj e U, but lj + N > I. 

It follows that 

d(f(z ),J f ) > r m+N V d(z l2 ,Jf) m > Dd(z ,J f ) m , (34) 

where D = r m+N r]C- 2m . 

Thus, in both case 1 and case 2 

d(f(z ),J f ) ^Dd(z ,Jf) m . 

It follows that Lemma H6] holds for S = D and K = max{KQ,m}. □ 

Theorem 47. There exists a TM which computes the coefficients of a poly- 
nomial p(n) such that for any z with d(z, Jf) > 2~ n one has 

f p{n) {z) £ V. 

Proof. Let zq G V\ and d(zo, Jf) > 2~ n . Let I be the maximal number, such 
that 

f j {z )eV 1 iorj = 0,l,...,l. 

Construct /j, 2$ as in Lemma I4TR Then properties 1) — 6) hold. Let s, I be as 
in Lemma 1451 and ii, . . . i s be the sequence of indexes for which 

Zi i v 1 (c\ u u. 

Then, by Lemma 1451 and properties 1) and 2) from Lemma 1461 

R > d{f{z), J f ) ^ C- 2 t s d{z, J f ) ^ C- 2 t s 2~ n . 

It follows that 

s < n\og t 2 + 2\og t C + \og t R. 
34 



Let % be an index such that z% G U 1 [C\ and U + No ^ I- Let 

fc^minOeN:/^^)^}. 
It follows from Lemma H2] that 

d( Zi , J f ) > 2- Kn , and 5 > d(f ki ( Zi ), J f ) > q k *d( Zl , J f ). 

Therefore, 

ki ^ K 2 n 

for some number K 2 > which can be obtained algorithmically. Now it is 
easy to see that among any Nq + K 2 Nn consecutive integer numbers between 
and I there exist i such that 

Zi £ U 1 (c\ U U. 

It follows that 

/ < {n log, 2 + 2 log, C + log, R) (N Q + K 2 Nn) . 
Since f^iV) m Vi m V m f(V), this finishes the proof. □ 

3.2 The algorithm. 

Similarly to Proposition 1331 one can prove the following result. 

Proposition 48. There is an algorithm computing two dyadic constants 
Ki,K 2 > such that for any z G V and any k G N if f k {z) G f(V) \ f~ l (V) 
then one has 

Ki K 2 

Now we are ready to describe the algorithm. The steps of the algorithm 
are analogous to the steps of the corresponding algorithm for a subhyperbolic 
map. Compute dyadic numbers s,e > such that 

e < wm{d{V 1} C \ f(V)), d{V } C \ V 1 )}, f(V) c U s (J f ). (35) 

Assume that we would like to verify that a dyadic point z is 2~ rt_1 close 
to Jf. Construct a dyadic set U 2 such that 

J f c U 2 m f-^V). 
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Then we can approximate the distance from a point z £ Ui to Jf by the 
distance form z to f/ 2 up to a constant factor. 

Now assume that z G Ui. Let p(n) be the polynomial from Proposition 
I4T1 Similarly to the case of a subhyperbolic map, consider the following 
subprogram: 
% : = 1 

while i ^ p(n + 1) do 

(1) Compute dyadic approximations 

Pi « /*(*) = f(r\z)) and d, » |£>/*(*)| = \Df~\z) ■ Df{f-\z))\ 

with precision min{2 -n-1 , e}. 

(2) Check the inclusion pi 6 Vj.: 

• if ^ G Vi, go to step (5); 

• if Pi ^ VI, proceed to step (3); 

(3) Check the inequality ^ fr22 n+1 + 1. If true, output and exit the 
subprogram, otherwise 

(4) output 1 and exit subprogram. 

(5) i ->■ i + 1 
end while 

(6) Output end exit, 
end 

The subprogram runs for at most p(n + 1) number of while-cycles each of 
which consist of a constant number of arithmetic operations with precision 
0(n) dyadic bits. The following proposition is proved in the same way as 
Proposition [3H 

Proposition 49. Let f(n, z) be the output of the subprogram. Then 

[ 1, if d(z,J f )>2- n -\ 

f(n, z)=l 0, if d(z, Jf) < K2' n -\ (36) 

y either or 1, otherwise, 

where K = J^) . , 
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4 Maps with parabolic periodic points. 



In this section we will sketch the proof of Theorem [TH Let / be a rational 
map such that flf does not contain either critical points or parabolic periodic 
points. We will assume that the periods and multipliers of the parabolic 
periodic points are given as a part of the non-uniform information. Replacing 
/ with some iteration of / if necessary, we may assume that the multiplier 
and the period of each of the parabolic periodic points is equal to 1. Then 
at each parabolic periodic point (in fact, a fixed point) p the map / can be 
written in the form 



where n p G N. We will assume that the numbers n p are also given as a part 
of the non-uniform information. Observe that some parabolic fixed points of 
/ may belong to the postcritical set of /. 

In this section we denote by CJf the set of critical points c G CJf such 
that the forward orbit O(z) does not contain any parabolic periodic points 
and by CJj the set of critical points c G Jf such that an iteration of c hits 

a parabolic fixed point. Let N (c), N , M, C be the same as in section [3J For 
numbers a, ft > 0, a direction v G [0, 2tt] and a point w G C introduce a 
sector 



Let pi, . . . ,p s be all parabolic fixed points of /. The algorithm will use 
the same non-uniform information (Nl, N2) which was described in sub- 
section [3J] and, in addition, 

N3. approximate positions cij of Pj with a dyadic precision a > and a 
number ft > such that for each j: 

1) pj is a unique fixed point of / in U2 a (cLj) (and thus pj can be approximated 
efficiently using Newton method); 

2) for each critical point c which lie in Jf one has 



n p +l 



+ 0(z n * +2 ), 



V^{ft, a) = {ze U a (w) : arg(z - w) G (u - ft, v + ft)}. 



(37) 



U 2a (a j )\{p j }nO(c) = 0; 
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Figure 6: Illustration to non-uniform information N3. 
belongs to an attracting Fatou petal at pj, 

W nJ f = M 

and the set 

u 2a ( Pj ) \ ((uv/) u { Pj }) 

belongs to the union of the repelling petals at pj. Observe that for a point 
z in the sector the distance d(z, Jf) up to a constant factor can be ap- 
proximated by | z — pj | . We will use the following result from [7] (see Lemma 
8). 

Lemma 50. Let g(z) = z + c n+ iz n+1 + c n+ 2Z n+2 + . . . be given as a power 
series with radius of convergence R > 0. There is an algorithm which given 
a point z with \z\ < 1/m < R computes I = [m n /C]-th iteration of z and 
the derivative Dg l (z) with precision 2~ s in time polynomial in s and logm. 
Here C is some dyadic constant which can be algorithmically constructed. 
The algorithm uses an oracle for the coefficients c n . 

Lemma [50] implies that we can compute " long" iterations of points close 
to parabolic periodic points efficiently. 

By the assumptions on the map / one can algorithmically construct a 
dyadic number e > such that the following is true. Let a point z be e- 
close to the Julia set. Assume that some iteration f k (z) belong to U e (pj). 
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Consider the neighborhood V = U e (f (z)). Let 



V = U{z),V u ...,V k = V 

be the pullbacks of V under / along the orbit of z, where U(z) is a neigh- 
borhood of z. Then the number of n for which V n contains a critical point 
c G Jf is bounded from above by a constant independent from k and z. In a 
similar fashion as Proposition [33] using ideas of Lemma HD] we can prove the 
following statement. 

Proposition 51. We can algorithmically construct dyadic numbers Ki, K2, e > 
such that for each z G U e {Jf) if k satisfies the following 

f n (z) G U e (Jf) for n = 1, . . . , k, and f\z) G U £ ( Vj ) 

for some parabolic fixed point pj then 

KW^Jf) ^ ^ K 2 d(f k (z),J f ) 



\Df*{z)\ ^ v ' /7 ^ \DP(z)\ ■ 

From the description of the dynamics near a parabolic point using Koebe 
Theorem [19] we can obtain the following result. 

Proposition 52. One can algorithmically construct numbers K3, L > such 
that for each z G U £ (pj) \ (UV^) and each n G N if 

f(z) G U E ( Pj ) foralli = 0,l,...,l = 2^ 
then one has \ f l {z) — Pj\ ^ 2 n \z — pj\ and 

Kl \f l (z)- Pj \ ^ d(f(z),J f ) ^ Ka \f\z)-Pj\ 



Pj\ ~ d(z,J f ) ^ ' \z-pj\ 

Now we briefly explain how to adopt the algorithm from Paragraph 13 . 21 to 
prove Theorem [TH Assume that we want to verify if a point z is 2 _n -close to 
Jf. Construct a sequence Zi = f li (z) in a way analogues to the construction 
from Lemma [46] (see f[32l ). We will define in a different way from f[32|) 
only if Zi G U e (pj) for some parabolic fixed point pj. In this case we do the 
following. 

1) If Zi G Vf for some attracting direction v then we stop. We can find the 
distance from Zi to Jf up to a constant factor. Using Proposition EH we can 



39 



estimate the distance from z to Jf up to a constant factor. 
2) If Zi £ U £ (pj) \ (UV^)) then consider the points f kr (zi), where k r = 2 [Lr] . 
Observe that by Lemma [501 we can find a 2~ n approximation of f hr (zi) in time 
polynomial in r and n. Let r be the minimal nonnegative integer number 
such that 

!/*-(*) -Pi | >e. 

Set = 

Observe that a direct analog of Theorem H7] is not true in a presence 
of parabolic periodic points, even if we assume that z does not belong to 
an attracting basin of a parabolic periodic point. For a point 2~ n close to 
a parabolic fixed point in a repelling petal it takes exponential time in n 
to escape an e- neighborhood of the parabolic fixed point. However, using 
Lemma EH and Proposition [52] we can prove the following. 

Proposition 53. There is an algorithm computing coefficients of a polyno- 
mial p{n) such that if 

d(z, Jf) > 2~ n and z { ^ UV? for each i = 0, 1, . . . ,p(n), 

then d(z p (n),Jf) > s. Moreover, we can compute 2~ n ~ 2 -approximations of 
z i — f ( z ) an d Df li (z) for i = 0, . . . ,p(n) in time polynomial in n. 

Now we briefly describe the algorithm computing the Julia set Jf in a 
polynomial time in n. Let z £ U e (Jf). Assume that we want to verify that z 
is 2~ n close to Jf. Without loss of generality we may assume that 2 1_n < e. 
Compute approximate values dj of Z{ and d{ of Df li (zi) with precision 2 -n ~ 2 , 
i = l,...,p(n). 

1) If Zi £ UVf for some i such that 

d(zj, Jf) ^ e/2 for j = 1, . . . ,i 

then we can find approximate distance from z% to Jf. By Proposition loTj we 
can find d(z, Jf) up to some constant. 

2) If d(zi, Jf) ^ e/2 for some i, then we can find the distance d(zi, Jf) up to 
a constant factor. Using Koebe Theorem f|T9|) we can find the distance from 
z to Jf up to a constant factor. 

3) If neither 1) nor 2) holds then by Proposition EH 

d(z, J f ) ^ 2- n . 
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In conclusion, let us mention an alternative numerical method to calculate 
iterations of points close to parabolic periodic orbits. This method has been 
implemented in recent literature. To illustrate the method we consider a map 
with a simple parabolic point at the origin: 

oo 

h(z) = z + z 2 + a k z k . 

k=3 

Let j(z) = —1/z and F(z) be the germ at infinity given by 

F(z)=jofoj(z) = -l/f(-l/z) = z + l + a[z), where a(z) = 0{z~ l ). 

Denote by T the unit shift: T(z) — z + 1. In [12] and [13] Ecalle has shown 
the following 

Theorem 54. The equation tyoF(z) = Toty(z) has a unique formal solution 
in terms of the series (generally, divergent) 

oo 

*(z) = P \ogz + Y,CkZ- k . (38) 

k=l 

The series ^> gives an asymptotic expansion for an attracting and a repelling 
Fatou coordinates of the map F , ^f a and ty r correspondingly. 

For survey on divergent series and asymptotic expansions we refer the 
reader to [20]. Theorem [541 means that the Fatou coordinates \l/ a and \EV near 
infinity can be approximated by finite sums of the series (155]) . To approximate 
n-th iteration of the map / near the origin one can use the following formula: 

f n ( z ) = J or 1 oro$o j, 

where ^ stands for either \l/ a or \l/ r depending on whether z belongs to an 
attracting or a repelling Fatou petal of /. We would like to emphasize that 
this is not a rigorous method since we do not know how many terms of the 
series (I38p to take to obtain the desired precision. However, empirically, the 
asymptotic expansion (138]) approximates the Fatou coordinates ^ a (z), *&b(z) 
with a very high precision. For instance, in [15] Lanford and Yampolsky 
used the series ( 155]) in their computational scheme for the fixed point of 
the parabolic renormalization operator. Our work in progress [11] give us a 
reason to hope that the described method can be made rigorous. 
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